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This is the Final Tcchnicil Report on NASA Gi'ent NSG H05 entitled 
'•Dyrtardo Inst.ibiHty of IXicts Ccnvoylni: Fluid," » rinp the ten-month 
period frrm October 197Ji to August 1975, awarded I!AbA I^glcy Research 
Center to Polytechnic Institute of Hew York, with Visiting Professor 
Yi-Yuan Yu as the Principal ^nvestirater. 

In this repot ♦ a ivry^:irj of the main results is presontrd first, 
rrose;.taticns ai;u puulicaticns oased oa the results arc ti;en itotnized. 

Three of these i are attached to this report ar Appendices. 

Staw'ary of Results 

There have been four separate, althourh closely related, research 
tasks cn ducts crn\cyit;p hdrh-eievd fluid. H'o main results are as follcrrs 

1. Crntilovorca Curved Ducts . The effect of curvature cn the flutter 
of cantilevered ducts har been examined i.n Rreat depth by conrtructinr root 
tocus dtpf’rcr*^ . Crr” rcheuri reru‘'ts hr.’"' hecr c'-t'ijred for a total of 
six ca-es, each inclueir.g the four ^otroct modes. Ar.onp other things, tlie 
results indicate that instaW ity of a curved luct can take place in the 
third node rr well a- in t)ie Sfccnd. Tn the care cf a str<al'’ht dvet. In- 
stability har b*?:C!'. ffurid previcus authers tt' -.ccur in the second mode 
only. 

2. Flexibly Sup]x.rted Ducts . To provide ari understanding of tiie 
d:*nar.ic in'- tr.bilLty ci.aractorlstics f real ducts which ire neither tixly 
cantilevered ncr aijr.^ly supfOi’tcd cr cL'ir.ptd, flexibly supported end con- 
ditions am err: id ;rf d. The rrlticai 'relccity as a fuixtion of the support 
flexibility har been sralyred. Corruter rmpra’^s for determining the 


criMc."'.! velccJty ar.:! for co \;tructir.f; rr.ot locwr dlai'nir.j wt ro rrcfar-rd for 
very I'.enori’ caroa of flcaciblo sup'-orta. For the special case of a straight 
duct with ono t^nd rigidly damped and the othor end flexibly aupj^ortedf 

conprchonsl VO nunorlcal results were obtained. The flexibility of the second 

. • 

ond was found to determine whether initial d^naanic instability woulu take 
place in the form of divcrrcr.ee or flutter. 

3. Finite Elenont. .•tnolyala . Hio dynrjaic instability analysis of arbi- 
trarily shaped ducts can bo carried out on tiie bacii of the finite elcsrent 
n.et3iod. Son-j oxdoraboiy work on tho finite element formulation of tho 
problen has been carried o\:t for the divergence and flutter of a otraight 
duct. Com'utcr proriran? have been prepared with the entire duct talren as 
one or two finite ele;r.onts. A llirlted amount of nmerical work has also 
been carried cut. 

Ij. IXict Systervj . For duct syctens cen osed of straight and circular 
aepr.cnts, another apprcach tc the djTianic instability problc.". is to consider 
the ducts as cenresite duct systtns and analyse them b;/ solving the appro- 
prl.aic Hiffercn*i'il r' *nt^cr.s and r/itchin- the beur.c;r.ry ccnditlona between 
adjacent segments. As an example, a ploiie syste:; censiscing of a circular 
segment witi: a straight one attached to each of its twc. ends was ccnsicered. 
Tro frc'!’ cnc” e^'-ntirn of the rroMcr. has beer, forr'ilatcd, but cenruter 
proprai-5 for the problem jcrnin yet to he ^Titter.. 

Presentations and Publications 

1. 7h*. Principal In’/cctigatf r presented the pJ^rer, "Application of 
Variational and GaJerkin P>qu?.tiens to Linear a:.d !!cnlinc?r Finite .ile.mcnt 
Analyc-'c," .at the Goncirss of t'-e Interrationa'' As+rcijautical Federa- 
tion in Ansterda.T., October 1971;. The work presented in the paper vas 
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initiated vh«*n he war Ir rjaldrnco at the VA5.A I-m-’lcy Research Center in 
the sunr.or j £ 1973. Tho procedure advocated in the paper has since been 
applied to the finite olonent fornmlation of the problem of dynamic insta- 
bility of straicht ducts, as described above under Task 3» 

2 , The Principal Investigator presented a Joint Seminar to the Depart- 
nont of Kechan5cal *!nf:ineerlnc and Department of Aerospace Engineering and 
Applied I'cchanlcs at Polytcchirtic Institute of Now York on October 31, 197U. 
The scrinar was entitled "Sore Recent N’Af /.-Supported Research.” In addition 
to a discussion of tho above application of variational and Galorkin equa- 
tions to tho finite element analysxs, the lecture also covered tho problem of 
4yna.nic instability of duct.-» convoying fluid. 

3. A presentation entitled "Some Recent NASA-Supported Research in 
Str.’Ctvral Kechr.:iies” was made by the Principal '^nvostifator at the llith 
Midvfestem Vccha-uics Corferenoo at Neman, 'klahcna, Karch 2U-26, 1975. 

’..■hilo Jimilar to tho seminar described in Item 2, the presentation included 
a froat deal riore results cn the d^.'nanic instability of ducts conve:,’lne 
fluid. An abr*~a3t of the presrntr tien arp erod in the ]^lur.o of Abstrac ts 
published by the Conference and is reprinted in this report in Appendix 1. 

L, A not: entitled '’Flutter ot' ■'’•ar.tilevcrcd Curved Ducts C'^nveying 
Fluid, " authored bp the '"rincippl Invest! tor, has been accepted for p'jbli- 
catinn in !'echa: :!cs Research Co;uTuiications . A preprint is included in 
this report in Ippcnd'.x 2. 

5. The Principal Investigator prepared the paper, "Effect of Curvature 
on the Flutter of Cantilevered Ducts Conveyinf’ Fluid," for presentation at 
the 11th Intcm?tior.al symposium on b'raco Tcchnolory and Science in Tokyo, 
July 1975. Cinoe the work has not been published in the U.S., permissicn to 
piesc-nt the paper was not t:,iant-d. A copy of it is attached as Appendix 3. 
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Sobs Recenc. NASA«Supporc.id Research In Structural Mechanics* 

By Y I -Yuan Yu 

Visiting Profccscr ol Mechanical Er.Rlneerlng. Polytechnic 
Institute, of . cw Yoric, and LrlaC In^julslicd Protessor of 
Aeronautical £r.^i.r.carlr.«i, Wicr.lta State University 


So::s recent research carried out by the author under the direct or indirect 
support of NASA is briefly reviewed. 

Tho first topic belnK investigated concerns the Calcrkin and variational 
forculation of tho finite clement analysis. As is known, the adv.int.ige of the 
Calorkin xec'aod lies with the fact chat its use does not depend upon the exist- 
ence of a variational principle. .Recently the mnehod has received auch atten- 
tion in connection wich fin'*'o eleaent formulation. This unfortunately can lead 
to S 0 .T .0 difficulties, as illustrated ay considering the problem of two- 
dir.ins ior.al elasticity. Vhen a variational principle dojs exist, such dlffl- 
culcica can be rexoved by resorting to the variational equation. 

The second topic deals with t'.e finite element analysis of nonlinear 
vibrations based upon the variational equation. The finite clement approxima- 
tion may be considered as a variational approximation with respect to the 
space coordinate. A second v.irlaclonal approx im.it Ion involves that with respect 
to time. The resulting nonlinear algebraic rcuatlons can then be solved by 
..n iteration procedure. To demonstrate the moched, vibrations of beams with 
largo deflection arc uiscussed. 

The last topic Is on the effect of high-velocity 'uid flow on the dynamic 
Instabllicv of ducts. Whe.i a certain critical velo< ' -f the fluid in a duct 

id reached, Instability of the- auct can develop in t. >. orm of divergence or 
flutter. Previous autnor.s rave treated mostly straig.iw ducts There appears 
to ■ . only or.j publication on curved duct fluid systems, and only divergence 
discussed in that publication. In the present work the flutter of cantl- 
li^vcrcd ducts with co.'.stant but nonzero curvature is investigated. Results 

th.'it tho effect of cur.Mture can be very severe. 

Rrieflv nentioned is the relation which the above topics of research bear 
•III. tho Space Hriitle ungine. 

*^The first tvo topics of r -Sv-rc.! described In this abstract were initiated 
when the author was in residence ac the Institute for Co.mputer Applications to 
Science a.-.d a-ngineering of .VASA Langley Research Center. .The third topic repre- 
sents work inltlitcd by the author when he w.is with Rocketdyne. The work is 
toing continued under the support ot .’.'ASA Langley .Research Center. 


Appendix 2 

(To be publishud in I'cchurdci; Ivescorch Ccrgrauicaticns ) 


FLUTfi.! OF CA’.’TILEVEIUiD CUHVED DUCTS COIIVST/IUG FLUID 


Yi-Yuan Yu 

Folybechnic Inr.tituto cf Kevr York, Brooklyn, Mow York 
V/icliita State University, V’ichita, Kansas 


Introduction 


VTien a certain critical velocity cf the fluid flew in an clastic duct is 
reached, instability of the duct con dcvelcp In the fom of diveri’onco or 
flutter, depend^ r.:; rainly or* th.e end conditions. Much has been done on the 
divergence and flutter of stmirht ducts; a comprehensive revie'.; u'as recently 
pd-ven by x'aido’:ssit ani Iscii Cll, Some pi’.bli cations on curved ducts have 
also appeared, bat noably dcalinn ulth diverper.ee. only Chon (2l further 
ur.alyr.'-d the oub-of-plan; flutter cf a cantilevered curved duct. In this 
work the in-plan-- flutter o.' a cantilevered cur*;od duct is investijatsd. 


Analysis 


iho duct has a lenfth L and radius a and contains a fluid flov; with a steady 
\-elocity 7. Fron Maniltcn's principle find the cor.plobe system ''f •• i- 
lincar equations of notion and associated toundarj' conditions, Ey *..71 ting 
frer. these the equations qcverring the irsLtial and fi-nal states and taking 
the difference, lir.sarised eo'aaticns governinp the psrturbed state are de- 
duced. 'Tnesc are further sir.pllfied by assuring inextcr.sional defemation 
and neglecting the tannentiai inoi ‘.ia. The assunpticn of inextension also 
climnates the effect of initial stresses. The sir.plificd equations are 
then used in the stability analysis of the duct-fluid cysten. As exact 
r.othods of solution are not feasible, the equatiens are solved by an itera- 
tive procedure cn the digital ccnputcr. 



Kumrlcal Rosulta 


Cfr'prohonuivo mmerical results have been obtained for the dir.ensionlcss cri- 


tical velocity ■ $/ I!/EI and frequency W ■ j/oI+nVil J1 L^, where M 

c c c c 

is tlic mass per unit length of the fluid flov;, n that of the duct, El the 
flexural rigidity of tho duct, and A the circular frequency# Those arc shown 
in Figs, 1 and 2 as functions of the mass ratio /3 ■ l!/lK+n) and tho curvature 
parameter 0 ^ “ L/a. Tho results for a straight duct for ^rtiich ■ 0 agree 
idth thoLo given bj' previous authors. Results for > 0 show that tho efiect 
of curvature can be veiy severe. The dimensionless ccmplex frequrncy eO has 
also been computed for a nvmbor of cases. Tne results snow that the third modo 
can become dynar.icallj' unstable first, in contrast to tho case of a straight 
cantilcv :rod duct for which tho second modo is known to become unstable first. 



FIG. 1 

Dir.cnsicnlcss Critical Velocity 



Dimocnsicnless Critical Frequency 
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a cimvATtmt on 

rUtTTia Of CAWIIiVLWD 31 JCT 5 CCHV^W', fUJTD* 
Tl-Tu«n TU** 


Abatnet 

ooily fluid flow In duet un Im> 1 b> dynamic InntablUty. Tbo 
rmnont study dMls with th« effaet of curvature on the fluttar of oantlle- 
verad d'icta. It was siotlvatad by the ••verely curved eonflfuratlon of the 
liquid fual and oxycan 'luota Inatalled on the apaca shuttle nain eniclne. 

1. Introduction 


When a certain critical valoelty of tha fluid In an ela.ftle duct Is reached, ln> 
itoblllty of the duct can daaalop In the fern of dlvcrpanoe or flutter, dependlnr > n 
Its end conditions. Ashlay and Havlland (Ref.l) wore tha first to intmatlfate the 
hvndlnx vlhrstlons of a 'laply supported duct. Housner (Ref. 2) re>exan!ned the pr-- 
blesi and showed that dlvcrKaneo can take place. Kls results were veriflod axperl- 
mentally by Dodds on<1 Runyan (Ref. 3). for cantilevered ducts c^nveylnr fluid Insta- 
bility can occui In the fnm of flutter. Ibis was first studied both analytically 
and experlmonta' ly by RenJaHln (Ref.lt) by considering a syste* of articulated rlftid 
ducts. Hla study wa'' later extended by Orefory and Paldouasla (Ref. 5) to a c''ntinu- 
ous el4ntle caatllivured duet. A ecaprehensl ve review on the dynamic instability 
of *tra: «tit ducta was recently Riven by Paldoussls and lasld (Kef. 6). 

Curved d’icta eonveylnit fluid have been InvoatlRated by Svetlltakll (Ref. 7), "nny, 
fartin, aod >iboy (Ref. 8), Chen (Refa.9, 10), and Hill ami Davis (Ref. 11). All '-f 
Vieae autivirs dl^Miiiaaed the divernence of curved ducta, and dvetiitskll and Hill and 
Davis further included the effect of Initial stresses. Only Chen (Ref. 10) also 
analyr.ed the out-of-plane flutter of a cantilevered duet. 

In this work the In-plane flutter of a cantilevered curved duct la Inveatipatel. 
Cosq:<l<de equatlrna IncludlnR tha effect of Initial stresses are first derived frov 
t'amllt'’n's principle. Tbese are rmt almplified. Rased on the simplified oquatirns 
the flutter analysis Is then carried out. ihsaerleal reeults are finally presented. 


?. Kqu .-i * jonn Of Motion and Boundary Conditions 

Ttic ctntlleverod circular duct under consideration is sh- wn in Fip.l, wt.are L ard 
a are the length and unlfoiei radius of curvature of the duct, reorectl voly, and V Is 
tho rtowlv voloeity of the fluid In the duct. The 
cooriMnites s, 0 am rnlated by s ■ aO. Applied to 
t'le t< tal Irnrtii of tho duct this gives I, • aw^,. 

For a duct with given L and a the total angle Oj, 
nay be consldarel t4> be the dlmenalonless curva- 
ture. 

To derlvo th<* eqiiatlcns of srtlon and tiie 
at stinted ticiin!:iry mndltlons f'T the duet-fluid 
stator, we uso M.viilton'a principle in tho form 

aO ( 1 ) 

The duct in taken to be the basic free body whose 
kinetic and strain energlee are T and U. The term 

4W then consists of the virtual work of all ex- Fig. 1 Curved Cart* lover 

*i?ma'< forces, which In this eaue are those exerted Ouet-Fluld Srtstem 

bv the fluid on the duct. The flnid la assumed to 



•T*i1s w>rk wa’ Initiated when the author was a consultant to Rocke’dyne. It wii 
.t'jh'»equentlv supi<>rted by NASA langley Research Center under NASA Orant ?!30 1105. 

'*Y!fltlng iTofessor rf Meehanleal .Jiginoering, Polytechnic Institute of New York, 
Rn ^klyn, New York, H.S.A., atd Distinguished Imfessor of Aeronautical r.ngineoring, 
WlchJta itate University, Wichita, Kansa.i, D.S.A. 




tv inr< r; ri'u: Iblc a<> ti>At nr intcmAl aitAiry !• MU'ciatvd wlUi It. 

i'nly ln>('Une Motlun 'T th« duet wlU b« oonnldaruJ, ‘Tie nrn*Mr> dlaplaewtnnt 
e<«t«'>n>mt. ' nt an arbitrary point in th« duet are thus thos« in Um radial and elrctaa* 
firentlal dlreeiluna. Pollowinr th« Wlaklar theory (H»r.l 7 ) im tajc« theaa in th* 

M, O* I/O JV ^ 2 ) 

wh- T'- u an>l v ar" Uie radial and elrr nforontlal ooMponanta ■ f the dlaplacaoont of a 
rclnt ' n the e''nt,roldal axla of tha •!. et, fron idilch t la waaured. A frlM dortotae 
llfrerontiatl''n with raapaet t<> a. 

kinetic oneriQf of tha duel la 


I 4idA 


w* >.r« ^ la the danally of tha duet Material, A the eroaa>aaettona1 area of tha duet, 
ari l a d-tt danrtoa dlf farantlatlon trlth raap^ to the tine t. Substitution frcai 
tqa.(2) ylaldJ . 

7 » j-/ 

uhori.' n la new tha na'’.i par unit lanrth of tha duet, IkUtlnp tha variation In tha ua> 
ual nan nr wa find . ri 

riTdt ^ 

avHI; ih«.‘ ftsju»Ad U|i And Up the only non-toro lirvi«r atrAln And rotation conponentA 

ATO 

t^a u*, 

and the eorrasrondlnr nrnllnaar atraln eonponent la 
The variation of tha strain anarr^ It then 

• -MjO Al% (M, ^ t - ^•‘Vj 


■dien* 


are the tension and bendln( 
traiiifoivtr. 1) Into 


t In tha duct, Dlffarantlati' n by rarta finally 


*-0 *•' ** 

Pi” fluid action on tho duet la due to the Inertia forces Induced by *he fluid 
accel or itlon. The conponents of the fluid acceleration in the rallal and eircun- 
farontlal dlraotlons have bean known (Ref, 11). If N la the nasa per <anli lenrth of 
tho fluid strean,the virtual work of the Inertia frrcaa on the duet la sijoply 

jW ■ - I''- 

~ j^M[ if i^)]4v4S (6) 

with IiQt.(?), ($), and (6) subsUtuU-d In ’laMllton'i principle In Ei].(l), wc find 



MV( 

fr'-r vtlcti thr> roilfirtx i^at«r of norCImar •qu.itlena of nutlcn an<l the aM'-clat-'d 
l<unu\ry a>ndttl''nr f»r tha oantllatMrad duet can ba mrlttan inaedlately. 

% Initial wd UTt>irt)ad Statae 


An liilttal utatc nf airaaa and daf' matlon of tha duct la rnaehad after a ctca<i/ 
ri'w la ratah'idicd In tha duet, Parturbatlm fr<jai U.e Initial atata will ba irtuoed 
by any fUrtla-r dlaturhanca of tha ayataai, uhleh nay bae<aie unatabla. Thtr final (lata 
C vemed i<y Kq>(7) ia thua equal to tha aun of the Initial and pcrturb>d atataa. 

With thtf mihKrlt ta n and a danotlnc thaaa atataa ua at^l write 

Ma«4«ai4«, vaVeft/* (6^ 

vnuatlcni rryiTnlnt* the initial atate can alan be written readily fn>r Ea.(7). ‘ 4 ?ua- 

tirra rnwmln> the lerturhad atata are than derived by aubatra.tlnr tha aquatl< nt 
frr the initial atata from thoaa for the final atate. The raaulta an- farther llnr- 
arlaed to yield 

uJi *1 d (imaj 

Uat; 

where and are fror. the aolutlon of the initial atate. 

The relatione betwnen Ma,Ra end Ua« ^a Forturbed atwta ate asriMu tu be 

ttic sane .If in tlio uaual static altuatlon: 

-6A2(a « ♦«*’«<; (Hi 

Nl'urc / is a dla»n8lt<nleB» constant of the Winkler theory defined by 

4A • -2A 

Kor a luet f r which a» r we have approiinately Z ■ I/Aa^, 1 bainr the nnnen<. of 
Im rtla <-f the cross aoetlon of the duet. 


ii. SlnpUflad bciuaijonL 


'r'nr vibration analysis of thin elreular linKr. Evrncen (Rnf.13) has ahown that 
■iccorocy In not. Inpalred by aeaunlnf Inertansl onal defcnantlcn (u^ * av^' > 0~' and 
l>y nerJeettns Ui" ttnnsnt.lal inertia (V* ■ 0). Inters:' tlnfly. the assunvtion of 
Iniortom. i'CvO ilfforwi' Icn will also ellninatc the effect of ttie Initial atreeser K, 
ar:d Nq. With Uio adoption of those assunptlons the second of Eqs.(9) beetiaes Iden- 
tically satisfied, and tho flr-'t of these reduces to 

(diSiL) (12) 

• * ^ * 

Tl"- v^a-tern in Kq.(lV') Is now asnuned to be neyliiclble In conparlson with U.e u,'« 
*.em praeedint’ It by adopting the often used practice for devcloplnr eqii.itlons of 
eyllnilrleal sheila (such as the well-known Donnell equations), which is to anrod- 
rute the rotation oopipontnt in 14. (li) ty lO- • -u*. With the subscript a further 
iropped we nay now write Eq.(l?) in tho following fornt 

(CtSiL) 113) 


lY.i- ass' elated tx iindary oonditlons for a eantllevared duct are now, 
vio; and (U), 

USU'aO u'*a*U*“a 0 


according to sqs. 

(IL) 


<S»LJ 



Tn «»v ->r»rtnul 'Jqr.nfi) t)w> br'tmriary ecndltlrn « Q at ■ • 0 ia !»>■ !■ nrrr nc< 1«<], 
an«l til- eoivliv *i K«/a • • 0 at • ■ I, la new Idantleal.y a«tlafltMt, 

Thf al«rliflr<t 'yatna of £^£.(13) and (Xi«) aiay alae b« d•rl^ed dlrartly fr*** 
l aaUU ti'a prlnelfla hr iaklnK 







Th«- af« V* rxirosalon f < r U har baon uaod by TlaMshanlib (Pmf.ll.). Tor a thlr rli • 

(7 • l/U') aiv! wltliovt fluid n<w (K • 0) l-^.(13) raduc-a to Uia ba»le aq-atl* u »*d 
by tvon^vn. For tan CMrvatura Eqa.(13) arul (lu) radvea to tha ntandari aq'iat: la 
of a *tra!fbt eontllenirad duet oontaylnr fluid, 
tlprii intn^elnr the dLranalonlaaa pararotarc 


ljqa.(13) and (11.) hoo<'ra 


i • */L, •) a y/u 

'"•/ft'"' 


Matt 


I0*]th (15) 


(T»U 

The •'tablllty annlynla to ha praaanta<l balo*» ic hatred on th» al.«rUflad a"uaflon!. 
Stability analyala baaad on the cnnplata aquatlona Mill ua earrlaJ out >aparataly. 


(V') 


’i. Stability Analyala 

For stability analyala eonaldar notion Klvan by tha etpraaelon 

>) * A (17) 

tdioro tt) la the dlronaionleaa fra<iuancy daflr-ad by M ■ 4 / {H * n)/£I J1 L* , fL brlnr 
ihu elrc«tlar franuoney. In faneral M la a eonplaz nvatbor. Hie ayatan wllJ ba 
.'itable or iinatahla dapandlni; upon uheth-'r tho lauflnary part of au la poaltlve or 
neratla^. For neiitral atablllty tho Inarlnary pari of U la taro. 

For a rrlvnn duet-fluid ayaten tha aiaar ratio 0 and dlneniilonloaa curvature are 
oo::lly enlculatad. Whan a value of the dlaionalohleaa flow valoelty v la further 
no1eei«o, an Infinite set of values tUn (» " 1, 2, can be found, alnoo the ays- 
ten posseeaos an Inflnlta nunber of darreea of fraadom. For v • 0, *Un real and 
arc thi dlmonoii nloaa natural freouanciev of the duct filled wl th a stationary fluid. 
Kn V Inereaseo, at>n will bacone eonplex, but tha ayster remains stable as lonr a.’ the 
Inarlnary part of A>n lo i«altlve. 1h:s la true until a cartalii critical vnlue v ■ 
vne In rearhed, when the lawfdnary part of (4',, bacoaiee terr araln, and a onndlilrn of 
neutral debility la arrived at in the nth lacK.e. Vfhan v > vne • the sysinn oaci'fieB 
iinotAble In Uic ntl. wda. Tha lowest amnnr nil Vno !• then the dlmanrlonl ess criti- 
cal flow velocity Vq, and tha eorrwspondlnr dlatcnslcnlasa critical frequency la 
•ien<itw»d ly which la neeasaarlly real, 

ITii* analysis nerlni with aubiitltutinr Eo.(17) Into (15). IM' ylalde the charac- 
terlcile aquation 

which is quartle. Asaune that the four roots are all distinct and denoted t>y 

dlBdlj 

Tho cfimlete solution is than ^ 

•1 • (19) 

J./ 
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Vih«n i| la nnw MbMitutad froi C(i.(I 9 ) tntr Um u.uiaury eenditl'na ( 16 ), w» rind 

2 A:» 0 . 1 ttfA: » 0 

fit * i** * 

i di-iti) A: « f 0 , i J-i A: e**j M 0 

jt| * * j$t * * ^ 

f I • trn-trtvial tiolutlan to «ilat, the dotonalnant of the ecofflelmta of nu^t 
*«nl3f . TMa river 
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«i 
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*1 
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^4 


0 <?') 


-.'|i>atl'ii (T>) 1« the oroii la froqueoey •quetlon, In Mhlch 4l j are funetlrni of the 
ilIrot"*'l -nless rrmiucncy Ul •!< iotomined by Eq. (Ifi). lilnee 4 la In feneral eoaqlci, 
) requlrau U>ni the reel and iMtrlnary parta of 4 vanlah te)>arataly. 

It la nni loa.iblo U> aolvo !qa.(l6^ and (?0) by direct awt.hoda, and lieratltra 
airtheih l/•r01h•r with the uac of dlrltal eoMputera have been reaorted U , The aqua* 
tl'du- Inir Ive the four [araaMtera ^ , 8;,, v, and ad, of which the flrat two (or any 
tve ) can hr rlvon ami the raaalnlng two are to be detemlne<l. In apflylnf the lte;a- 
• on rr ce<iure when $ and 8j, are riven, trial valuea of v and ad ar» fl^^t aas'jted. 
Kqiiatlon (Ih) la then aolved nuatarleally for *1 j, but true reau''ta of A, are 
o-ta1ne<l only wtien the real and ijaarlnary part a of 4 In hq.(?0) botli vanish. Thla 
It' aehlrved ny Itaratlni v and |8 on the ooaqniter. 


. . 1 , -rl cal Keaulta 

'cnrrrhriipi VO ntMTleal results hrve t'sen obtained for v^ and (d,. for the full 
f'lrro < f 0 i$t 1 ai'd for Bf ■ 0, 1, 6, aa ahown In Firs* ? ar.d ), The- 

results for a atralidit duet for whleh 8j. ■ 0 arree with thus# r1 veil bv | revloua 
authors, Reaultk for wj, > 0 ah' w that the effect of ourvature can i>e ver- sever*. 

In renrrit • Ve Increasea with Increealnr 8^ In the ranr* of larr* fi Oieavler fluid), 
(bi the ottier liantl. It decreaaes with Inereaflnii 8i, In the ranre of aina)] fi (llfhter 
Puld), isit there saesii. to be a lower bound near 8x, ■ 5, 6. 1h* result for We for 

the st.rairht duet always prevldas a loiror bound to the reaulta for curved duets. 

The value »f We Ineren'ea ■onotcnleally with 6v, until 8], ■ 1., tieyrinc' whleh the trrnd 
i'ccoaw'8 rntlHir irrerular. 




Klf. ? Idnenpl inlass Iritloal Velocity 
Ve la Function of Hasa Ratio 0 
an t lurvaturo KarsMater Bj, 


FiR. 3 Dljsenslonloas Critical Frequency 
adc aa Function of Haas Ratio 0 
anrl Curvature Paraan-ter 8[, 


V 


\ 



The 'ilMMionlt*. c’Hi'lm fr«qu*nc]r u) iuo i«an c«Ni«t«d fur a nuwtKr of eaeoa. 

Ttio r^ulta for tho lowaat four nodar in tho ease of it ■ 0.?$ and ■ 2 are ahowr. 

In flR. bt Mtilch la Mietitlally a root locua dlaftran. Aa nontlonod Mfora, when 
V * 0, M 11 real aii'l the branehaa InU-raect the horl rental axla at polnta correa- 
petdlnr to tho nttural fraquanelaa of a duct with a itatl’-naiy fluid. Ai v Incroasea 
W bocoMaa cofifiiax, and all four MMSa.i are (landed at flrvi. Aa v an rcaehea 6.5t 
for the third node t>eoonaa ra«. attain and la equal to affro-xlnately 3^ (U.eae valuer 
of V an* ta belnr wirlflable by the reaulta in fir*. 2 and 3). Ihia la tha dlner- 
alenlr.aa erltical velocity v.. When It la axceeda-i, tha duet baoona: unatabla. It 
In Ir’ereatlnc to .wta that the third node can become dynaatlcally unaiabic flrct. In 
contrac t to the cane of a atraifiht cantilevered duct for which tha acc< nd node haa 
been known t<< bceonc inatablc firrt. 



Mr. b Root lioeun Dlarraai (0 ■ 0.35, Sj, ■ 2) 
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